Engineering arbitrary pure and mixed quantum states 
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This work addresses a fundamental problem of controllability of open quantum systems, meaning 
the ability to steer arbitrary initial system density matrix into any final density matrix. We show 
that under certain general conditions open quantum systems are completely controllable and propose 
the first, to the best of our knowledge, deterministic method for a laboratory realization of such 
controllability which allows for a practical engineering of arbitrary pure and mixed quantum states. 
The method exploits manipulation by the system with a laser field and a tailored nonequilibrium and 
time-dependent state of the surrounding environment. As a physical example of the environment 
we consider incoherent light, where control is its spectral density. The method has two specifically 
important properties: it realizes the strongest possible degree of quantum state control — complete 
density matrix controllability which is the ability to steer arbitrary pure and mixed initial states 
into any desired pure or mixed final state, and is "all-to-one", i.e. such that each particular control 
can transfer simultaneously all initial system states into one target state. 



Controlled nianipulation by atoms and molecules us- 
ing external controls is an active field of modern re- 
search with applications ranging from selective creation 
of atomic or molecular excitations out to control of chem- 
ical reactions or design of nanoscale systems with desired 
properties. The external control may be either coherent 
(e.g., a tailored laser pulse or incoherent (e.g., a spe- 
cially adjusted or engineered environment 0, Q or quan- 
tum measurements commonly used with and sometimes 
without feedback 

A challenging topic in quantum control is to pro- 
vide practical methods for engineering arbitrary quantum 
states [H. The interest to this topic is driven by funda- 
mental connections to quantum physics as well as by po- 
tential applications to quantum state measurement [Ej 
and quantum computiiig with mixed states and non- 
unitary quantum gates [6] . Various recipes for engineer- 
ing arbitrary quantum states of light were proposed 
For matter, engineering arbitrary open system's quan- 
tum dynamics with coherent control, quantum measure- 
ments and feedback was shown to be achievable ^ . How- 
ever, the problem of deterministic engineering of arbi- 
trary quantum states of matter has generally been re- 
mained unsolved. 

This works proposes the first, to the best of the au- 
thor's knowledge, deterministic method for engineering 
arbitrary pure and mixed density matrices for a wide 
class of quantum systems. The method uses a combina- 
tion of incoherent control by engineering the state of the 
environment (for which we consider appropriately filtered 
incoherent radiation) on the time scale of several orders 
of magnitude of the characteristic system relaxation time 
Tici followed by fast (e.g., femto-second) coherent laser 
control to produce arbitrary pure and mixed states for a 
wide class of quantum systems. The method is determin- 
istic in the sense that it does not use real-time feedback 
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and can be applied to an ensemble of systems without the 
need for an individual addressing of each system. Impor- 
tant is that the suggested scheme requires the ability to 
manipulate by both Hamiltonian and non-Hamiltonian 
aspects of the dynamics; a fundamental result of Altafini 
shows that varying only the Hamiltonian is not sufficient 
to produce arbitrary states of a quantum system [Ifl] . 

Engineered environments were suggested for improv- 
ing quantum computation and quantum state engineer- 
ing making robust quantum memories (l2| . prepar- 
ing many- body states and non-equilibrium quantum 
phases |13| , inducing multiparticle entanglement dynam- 
ics This work exploits incoherent control by engi- 
neered environment [2] in a combination with coherent 
control for producing arbitrary pure and mixed density 
matrices. While incoherent processes were used in var- 
ious circumstances, as for example in cold molecule re- 
search, to prepare a pure state needed for a full control 
over the system's pure states [Tsj . their use in the pro- 
posed method serves for a more general goal of engineer- 
ing arbitrary pure and mixed quantum states. 

The method has two special properties. First, it im- 
plements complete density matrix controllability — the 
strongest possible degree of state control for quantum 
systems meaning the ability to prepare in a controllable 
way any density matrix starting from any initial state. 
Second, the produced controls are all-to-one — any such 
control c* transfers all pure and mixed initial states into 
the same final state and thus can be optimal simultane- 
ously for all initial states [l^ . This property has no ana- 
log for purely coherent control of closed quantum systems 
with unitary evolution, where different initial states in 
general require different optimal controls. While an ab- 
stract theoretical construction of all-to-one controls was 
provided jT6| . the problem of their physical realization 
has remained open. The suggested method provides a 
solution for such a physical realization. 
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I. COHERENT AND INCOHERENT 
CONTROLS 

The dynamics of a controlled n-level quantum system 
isolated from the environment is described by density 
matrix pt satisfying the equation 

pt = -i[H^) + u{t)V, pt] , pt=o = Pi- (1) 

Here Hq = ^^^'^ system Hamiltonian 

(with eigenvalues ei < £2 < ■ • • < a-nd eigenvectors 
\i)) and V is the interaction Hamiltonian describing cou- 
pling of the system to the control field u{t) (e.g., a shaped 
laser pulse)l; commonly V = —p, where p is the dipole 
moment of the system. The evolution is unitary, pt = 
UtpoUj, where the unitary evolution operator Ut satis- 
fies the Schrodinger equation Ut = —i[Ho + u{t)V]Ut- 
The unitary nature of the evolution induced by the field 
u(t) implies preservation of coherence in the system such 
that for example pure states will always remain pure; the 
corresponding control is called coherent. 

If the system interacts with the environment, then its 
evolution in the absence of coherent control {u{t) = 0) is 
described by the master equation for the reduced density 
matrix (itI fisj 

Pt = -i[Ho + Hcs,Pt] + C{pt), pt=o = Pi (2) 

where the superoperator C and the effective Hamiltonian 
HcS describe the influence of the environment. The ef- 
fective Hamiltonian H^s represents spectral broadening 
of the system energy levels and typically commutes with 
Hq. For a Markovian environment the superoperator C 
has the form C{p) = J^C^^iP^l ~ L\Lip — pL\Li) with 

some matrices Li [l^ ■ The explicit form of these ma- 
trices is determined by the state of the environment and 
by the details of the microscopic interaction between the 
system and the environment. 

The matrices Li are usually considered as fixed and 
having deleterious effect on the ability to control the sys- 
tem — open quantum systems subject to the Markovian 
evolution ^ with constant C are uncontrollable 
However, the assumption of constant C is too restric- 
tive, since L can be manipulated by adjusting the state 
of the environment through its temperature, pressure, or 
more generally through its distribution function (spec- 
tral density). Control through adjusting the state of the 
environment in general does not preserve quantum co- 
herence in the controlled system and for this reason is 
called incoherent [1]. 

This work considers incoherent radiation as an exam- 
ple of a Markovian control environment and the scheme is 
analyzed below for this case. Other environments, either 
Markovian or non-Markovian, can also be used as inco- 
herent controls; however, the ability to use a particular 
environment for engineering arbitrary quantum states us- 
ing the proposed scheme requires a separate analysis in 
each case. The state of the environment formed by in- 
coherent photons is characterized at a time moment t by 



spectral density nk(t) of photons with momentum k; in 
general, spectral density of photons with polarization a 
can also be exploited. For the purpose of this work the di- 
rectional dependence of spectral density is not necessary 
and its dependence only in the photon energy cj = |k| is 
used; is assumed to be constant over the frequency 
range of significant absorption and emission for each sys- 
tem transition frequency uJij = £j — Si. The spectral 
density (t) can be experimentally manipulated for ex- 
ample by filtering. 

The evolution of the system in the environment formed 
by incoherent radiation with spectral density is de- 
scribed by the master equation ^ with superoperator 
C = Cn^ of the form (MIHI 

(3) 

where Aij > are the Einstein coefficients for sponta- 
neous emission, ujij = Si — £j are the system transition 
frequencies, Qij = is the transition operator for the 
\j) \i) transition, and Lq{p) = 2QpQ^ — Q^Qp—pQ^Q. 
The intensities u{t) and n^j of the coherent and incoher- 
ent light are the coherent and incoherent controls, re- 
spectively. Radiative energy density per unit angular 
frequency interval is p^j = fujjn^{uj'^ / 1:'^ c^) . 



II. ENGINEERING ARBITRARY QUANTUM 
STATES 

Let Pi be any (mixed or pure) initial state of the sys- 
tem and pf — ^Pi\4'i){4>i\ be an arbitrary (mixed or 
pure) target (final) state. Without loss of generality we 
assume pi > P2 > • ■ • > • The system is assumed to be 
generic in the sense that all its transition frequencies are 
different (in particular, its spectrum is non-degenerate) 
and all Einstein coefficients Aij are non-zero. The goal 
is to find a combination of coherent and incoherent fields 
transforming all pi into pf. 

The scheme to steer pi into pt consists of two stages. 
In the first stage, the system evolves on the time scale 
of several orders of magnitude of the relaxation time, 
t ss OTi-ci (where a can be chosen in the range 2 10 
depending on the required degree of accuracy) under the 
action of a suitable optimal incoherent control into 
the state pf = '^Pi\i){i\ diagonal in the basis of and 
having the same spectrum as the final state pf . The state 
Pf has the same purity as pt; pt is mixed if pf is mixed 
and Pt is pure if pf is pure. In the second stage, the 
system evolves on the fast (e.g. femto-second) time scale 
under the action of a suitable coherent laser control which 
rotates the basis of Hq to match the basis of pf . 

The first stage exploits incoherent control with any n*^ 
such that n*^^. = Pj/(j>i — Pj) to prepare the system in 
the state pf. Coherent control is switched off (u(t) = 0) 
during this stage and the system dynamics is described 
by the master equation ^ which for off-diagonal matrix 
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elements pin = ('IpI^^) of the density matrix takes the 
form 

Pin = -{iain + r/„)p;„, r;„ = ^{Wji + Wjn) > 

i 

where a/„ = {l\Ho + Hcs\l) - {n\Ho + Hes\n) and Wy = 
Aij {2ni^ ■■+!). Off-diagonal elements decay exponentially 
for r;„ > 0, pin ~ exp(— r;„f). Diagonal elements satisfy 
the Pauli master equation 

Pnn = y^S^njPjj ^ WjnPnn) 
3 

For generic systems, all rii^^. can be independently ad- 
justed, and the detailed balance condition WnjPjj ~ 
WjnPnn implies n*.^. — Pj/{pi — Pj)- The master equa- 
tion for spectral density such that n*^.. — pjj (pi —pj) 
has pf as the stationary state and exponentially fast 
drives the system to pf [2^. The case when some 
Pi = Pj formally requires infinite density at the cor- 
responding transition (e.g. infinite temperature of the 
environment is required for creating equally populated 
state Pf = il), but for practical applications a reason- 
able degree of accuracy allowing for a finite n* . . is always 
sufficient. The first stage when necessary can i^e divided 
into two parts, diagonalization of the density matrix with 
any producing positive and sufficiently large Tin, fol- 
lowed by the evolution produced by any n* such that 
= PoKPi - Pj) to produce pi. 
The second stage implements coherent dynamics with 
unitary evolution transforming the basis into {|(/>i)} 
and hence steering pf into pf. For successful realization 
of this stage we assume that any unitary evolution op- 
erator of the system can be produced with available co- 
herent controls when >C = on a time scale sufficiently 
shorter than Trei, i.e. that the system is unitary control- 
lable when decoherence effects are negligible. Incoherent 
control is switched off during this stage by setting = Q 
and the dynamics is well approximated by the unitary 
evolution ([1]). 

Necessary and sufficient conditions for unitary control- 
lability were obtained by Ramakrishna et al. [23j. Suf- 
ficient conditions for complete controllability of n-level 
quantum systems subject to a single control pulse that 
addresses multiple allowed transitions concurrently were 
established [24]. Ramakrishna et al. showed that a nec- 
essary and sufficient condition for unitary controllability 
of a quantum system with Hamiltonian H = Hq +u{t)V 
is that the Lie algebra generated by the operators Hq 
and V has dimension n^. We assume that the system 
satisfies this condition when £ = and that any U can 
be produced on a time scale sufficiently shorter than the 
decoherence time scale, e.g. by using time-optimal con- 
trol dl]. Under these assumptions a coherent field u*{t) 
exists which produces a unitary operator U transforming 
the basis {\i)} into {|(/>i)} and therefore steering pi into 
Pf. This field implements the second stage of the scheme 
to finish evolving pi into pf. The second stage can be 



implemented also on the long time scale using dynamical 
decoupling [2^ to effectively decouple the system from 
the environment. 

There exist other methods to prepare arbitrary quan- 
tum states. An example is the method of engineering 
arbitrary Kraus maps proposed by Lloyd and Viola Q 
which in particular can be used to produce arbitrary 
quantum states. Another option is to cool an ensemble 
of systems to a pure state |0) and then apply randomly 
and independently to each system a control producing a 
unitary operator Ui transforming into \4>i). Imple- 
menting each Ui with probability pi will effectively pro- 
duce the systems in the state pf = J2Pi\4'i){4>i\- Both 
these methods require independent addressing of indi- 
vidual systems in the ensemble that in general is hard 
to realize. The proposed method is free of this short- 
coming and can be applied to an ensemble without the 
need to independently address each system. This ad- 
vantage should significantly simplify practical realization 
when individual addressing of quantum systems is hard 
to realize. We remark that the first control stage re- 
quires individual addressing of each transition frequency. 
If the system is non-generic, in particular, is degenerate 
then the method may work if the degenerate states can 
be selectively addressed by using polarization of the in- 
coherent radiation. For degenerate systems (harmonic 
oscillator) another scheme was proposed Q . 

Determining the correct form of the master equation 
for systems with time dependent Hamiltonians in general 
is a non-trivial problem [23, H^. However, the proposed 
scheme does not use the master equation with time de- 
pendent Hamiltonian u(t)V and therefore the problem 
of choosing the correct form of the master equation with 
time dependent Hamiltonian is not relevant for this work. 
The scheme in its first step exploits incoherent control 
when the coherent control is switched off and therefore 
the system Hamiltonian is time independent. In this case 
the system dynamics under Markovian approximation is 
governed by a master equation of the form ^ . The sec- 
ond step exploits fast coherent control on a short time 
scale when incoherent control is switched off and the de- 
coherence effects are negligible. In this case £ ~ and 
the dynamics is well approximated by ([Ij. 

If the quantum system is sufficiently simple and all its 
relevant parameters are known, then methods of optimal 
control theory (OCT) can be used to find optimal fields 
n* and u*{t) as shown on Fig.[Tl subplot (a). If the sys- 
tem is complex and (or) some of its relevant parameters 
are unknown then various adaptive learning algorithms 
can be used to implement the proposed scheme for engi- 
neering arbitrary quantum states (Fig.[l] subplot (b)). In 
this case, the creation of a desired diagonal mixed state 
Pf during the first stage of the control scheme can be 
realized with learning algorithms The second stage 
can be implemented using learning algorithms for coher- 
ent control pgj . Learning algorithms were shown to be 
efficient for both cases and therefore they will be efficient 
when used in a successive combination as required in the 
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FIG. 1. (Color online). Implementing complete density ma- 
trix controllability with (a) optimal control theory (OCT) and 
(b) learning control. In (a), methods of OCT and a theoret- 
ical model of the system dynamics are used to find optimal 
incoherent and coherent controls n\ and u*(t). In (b), the 
objective J — \\pT — Pf||, where px is the output state, is 
measured (block with letter "M ") at the end of each iteration 
and the result is used to design a new control, ideally with 
better performance. The cycle is repeated until obtaining a 
satisfactory output. 



proposed scheme. 



III. COMPLETE DENSITY MATRIX 
CONTROLLABILITY 

Depending on a particular problem, various notions of 
controllability for quantum systems are used ^Ifli, il6|, i23i , 
[30l - [33 |. including unitary controllability, pure state, den- 
sity matrix and observable controllability [s^l, and com- 
plete density matrix controllability [l6l |. Unitary con- 
trollability means the ability to produce with available 
controls any unitary evolution operator U. Density ma- 
trix controllability means the ability to transfer one into 
another arbitrary density matrices with the same spec- 
trum (i.e., kinematically equivalent density matrices); a 
particular case is pure state controllability as the ability 
to transfer one into another arbitrary pure states. Com- 
plete density matrix controllability means the ability to 
steer any initial (pure or mixed) density matrix pi into 
any (pure or mixed) final density matrix pt, irrespective 
of their relative spectra [T6j . This notion is different from 
density matrix controllability (soj . where only kinemati- 
cally equivalent density matrices are required to be acces- 



sible one from another, and is the strongest among all de- 
grees of state controllability for quantum systems; com- 
plete density matrix controllability of a quantum system 
implies in particular its pure state and density matrix 
controllability. The suggested scheme allows for transfer- 
ring one into another arbitrary density matrices thereby 
approximately realizing complete density matrix control- 
lability of quantum systems — the strongest possible de- 
gree of their state control. 



IV. ALL-TO-ONE CONTROLS 

An all-to-one control c, is a control which steers all ini- 
tial states into one final state. Such controls can be opti- 
mal simultaneously for all initial states [l^ and their im- 
portance is motivated by the following. Let J(c) = J{pt) 
be an arbitrary control objective determined by the sys- 
tem density matrix px evolving from the initial state 
Pq to the final time T under the action of the control 
c (e.g., J = Ti[ptO] for some Hermitian observable O 
or J = IIpt ~ Pf ID- In general, optimal controls — con- 
trols minimizing the objective — are different for differ- 
ent initial states. However, if c* is an all-to-one control 
steering all states into a state pf minimizing J, then c* 
will optimize the objective for any initial system state. 
While an abstract theoretical construction of all-to-one 
controls was provided in terms of special Kraus maps 
whose definition is also provided below, their physical re- 
alizations has remained as an open problem. The pro- 
posed control scheme provides a physical realization of 
all-to-one controls and the corresponding Kraus maps for 
any pure or mixed state pf. The all-to-one property is 
achieved during the first stage, where n* produces the 
same density matrix pf independently of the initial state. 

The density matrix representing the state of an n-level 
quantum system is a positive, unit-trace n x n matrix. 
We denote by = ^nxn -j-j^g all n x n complex 

matrices, and by Vn := {p G Mn \ P — P > 0, Tr(/o) = 
1} the set of all density matrices. A map $ : — 
is positive if ^(p) > for any p > in A4n- A map $ : 
Mn — > Mn is completely positive (CP) if for any I e N 
the map $ (g) I; : Mn ^ Mi Mn <8) Mi is positive (1/ is 
the identity map in Mi). A CP map is trace preserving 
if Tr[$(p)] ~ Tr(p) for any p e Mn- Completely positive 
trace preserving maps are referred to as Kraus maps; 
they represent the reduced dynamics of quantum systems 
initially uncorrelated with the environment [isl [33| . 

Any Kraus map $ can be represented using 
the operator-sum representation (OSR) as $(/?) = 
'^^^i KipXj , where {Ki} is a set of complex n x n ma- 
trices satisfying the condition X]i=i ^l^i — to ensure 
trace preservation ^18, 33]. The OSR is not unique: any 
Kraus map $ can be represented using infinitely many 
sets of Kraus operators. 

The all-to-one Kraus map for a given final state pf = 
'^iPi\'t>i){4>i\ is defined as a Kraus map '^p^ steering all 



initial states into pf, i.e., such that '^pf{p) 



Pi 



for all 
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P € T^n [16|. If pf is a density matrix maximizing a 
given objective J and c* is a control producing the map 
<i>Pf , then this control will be simultaneously optimal for 
all initial states, i.e., the same c* will maximize the ob- 
jective for any initial system state. An OSR for a uni- 
versally optimal Kraus map $pj. can be constructed by 
using Kij = ^/Pi\4>^){4>3 \ ^ Kraus operators. Indeed, 

for any p, Ei^i ^^JPK^J = Pf ■ 

All-to-one Kraus maps were constructed theoretically 
in [l^. The two stage control scheme described in the 
present paper provides an approximate physical realiza- 
tion of all-to-one Kraus maps for all pf (therefore any 
all-to-one Kraus map can be produced using this scheme) 
for generic systems unitary controllable on the fast with 
respect to TicI time scale. 



V. EXAMPLE: CALCIUM ATOM 

The scheme is illustrated below with an example of a 
two-level atom whose all relevant parameters are known 
and it is easy to analytically understand and visualize 
the controlled dynamics. We consider calcium upper 
and lower levels 4^P and 4^S as two states |1) and |0) 
of the two-level system. For this system the transition 
frequency is UJ21 — 4.5 x 10^^ rad/s, the radiative life- 
time t2i = 4.5 ns, the Einstein coefficient A21 = l/t2i ~ 
2.2 X 10^ s"\ and the dipole moment p.12 = 2.4 x 10"^^ 
C-m [3^ . The method works equally well for any ini- 
tial and target states, we take for the sake of definiteness 
Po = |0)(0| andpf = i|0)(0| + ||l)(l|. 

The system is generic, all its relevant parameters are 
known and we can analytically find optimal controls. The 
goal of the first (incoherent) stage is to prepare the mix- 
ture Pf = ||0)(0| -I- This goal is realized by ap- 
plying to the system incoherent radiation with spectral 
density satisfying rit^^i — 1/2 during time T of several 
magnitudes of decoherence time ^21 ; we choose T — 50 ns. 
The goal of the second (coherent) stage is to rotate the 
state produced at the end of the first stage to transform 
it into Pf . This goal can be realized for example by apply- 
ing a resonant 7r-pulse E{t) = Ecos{oji2t). Other meth- 
ods of coherent control which produce the same unitary 
transformation can be used as well. The electric field 
amplitude which makes the Rabi frequency equal to the 
radiative decay rate is i5 « 10'^ V/m [3J]. The pro- 
posed scheme requires the duration of the second stage 
to be significantly shorter that the decay time. This can 
be satisfied by choosing E ^ 10^ V/m. We take reso- 
nant electromagnetic field E{t) — E cos{uJi2t) of ampli- 
tude E = 10^ V/m acting on the system during the time 
interval T{ ~ T = 1310 fs. The Rabi frequency for the 
field of such amplitude is I^r « 1/2320 fs~^, thus the field 
acts as a 7r-pulse transforming the state pf into pf. The 
results of the numerical simulation are shown on Fig. (3] 
The time interval 1310 fs is much less than ^21 and de- 
coherence effects are negligible during the second stage. 
For illustrative purpose of better visualization of the tra- 



jectory during the second stage we provide in Fig. [3] sim- 
ulation results for E = 10^ V/m and Tf - T ^ 13.1 fs. 
The method can be applied for producing arbitrary pure 
or mixed target density matrices from any pure or mixed 
initial state thereby implying complete density matrix 
controllability of this system. 




Time 



FIG. 2. (Color online). Application of the control scheme 
for engineering a desired mixed state of two Ca energy lev- 
els. Subplot (a) shows the behavior of the objective vs. time 
during the two stages of the control. The timescale is not 
uniform, the duration of the first stage is T = 50 ns and of 
the second stage is Tf — T = 1310 fs. Solid line is for \\pt — pi\\ 
and dotted line for \\pt — Pf ||; the latter quantity almost com- 
pletely vanishes during the first stage of the control. Subplot 
(b) shows the corresponding evolution of the Bloch vector. 
The initial system state corresponds to the south pole, the 
vertical blue line represents incoherent evolution during the 
first stage and the green line shows coherent rotations during 
the second stage under the action of resonant field of ampli- 
tude E = lO' V/m. The evolved state approaches the target 
state indicated by the red marker. 




FIG. 3. (Color online). The only difference from the descrip- 
tion of Fig. [2]is that E = lO'-* V/m and Ti~T= 13.1 fs. 
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